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We show that a quantum-limited phase-preserving amplifier can act as a which-path information 
eraser when followed by heterodyne detection. This Team splitter with gain’ implements a continuous 
joint measurement on the signal sources. As an application, we propose heralded concurrent remote 
entanglement generation between two qubits coupled dispersively to separate cavities. Dissimilar 
qubit-cavity pairs can be made indistinguishable by simple engineering of the cavity driving helds 
providing further experimental flexibility and the prospect for scalability. Additionally, we hnd an 
analytic solution for the stochastic master equation—a quantum filter—yielding a thorough physical 
understanding of the nonlinear measurement process leading to an entangled state of the qubits. 

We determine the concurrence of the entangled states and analyze its dependence on losses and 
measurement inefficiencies. 


I. INTRODUCTION 

Spatially separated objects can be entangled by the 
measurement backaction [1-3] of a joint measurement. In 
the process of concurrent remote entanglement generation, 
two remote stationary qubits are first entangled with sep¬ 
arate flying qubits and which-path information is erased 
from them by interference effects [4, 5]. A subsequent 
measurement of the flying qubits can then implement a 
joint measurement with a backaction that projects the 
stationary sources to an entangled state, demonstrated 
in various atomic and solid state systems [6-12] through 
coincidence detection of photons [13]. In this concurrent 
scheme, entanglement generation occurs purely by mea¬ 
surement backaction; the entangled qubits exchange no 
information, not even unidirectionally. There need be no 
causal connection between the qubits. This represents an 
important conceptual difference from consecutive remote 
entangling configurations [3, 14] with unidirectional ex¬ 
change of information due to both qubits seeing the same 
photon field. 

Generating entanglement is a necessity for quantum 
communication, cryptography and computation [15-17]. 
The concurrent configuration promotes scalability and 
modularity of a quantum network, allowing entangling 
operations between arbitrary nodes through routing inde¬ 
pendently generated parallel signals to a quantum eraser. 
Compared to consecutive configurations, high entangle¬ 
ment fidelity is harder to achieve since the single-qubit 
information is more exposed to losses. However, the 
concurrent method provides better on-off ratio for the 
effective entanglement since, with use of directional ele¬ 
ments, no parasitic signal could in principle propagate 
from one qubit to the other. 

We show that a quantum-limited phase-preserving am¬ 
plifier can be used as a quantum eraser for the which-path 
information for concurrent microwave signals (Sec. II). 
When followed by detection of both quadratures of the 
amplified output signal, this novel ‘beam splitter with 
gain’ can implement continuous joint measurement on 
the remote signal sources. The quantum eraser config¬ 
uration is general for systems involving continuous vari¬ 


ables [18, 19]. For concreteness, in Sec. Ill we propose 
and analyze remote entangling for superconducting qubits 
coupled to traveling continuous microwave signals (see 
also Ref. [20]). Analogously to the spatiotemporal mode 
shapes, the single-qubit information from dissimilar cavity- 
qubit pairs (unequal dispersive coupling or decay rates) 
is carried by the unequal temporal measurement ampli¬ 
tudes. However, because the measurement amplitudes 
depend on the cavity dynamics, they can be made indistin¬ 
guishable through simple engineering of the cavity driving 
amplitudes, reinforcing the scalability and experimental 
flexibility of the entangling scheme. In Sec. IV we derive 
an analytic solution for the qubits’ stochastic measure¬ 
ment dynamics and analyze the fidelity of the resulting 
entanglement with realistic estimates before concluding 
in Sec. V. 


II. PHASE-PRESERVING AMPLIFIER AS 
A WHICH-PATH INFORMATION ERASER 

Quantum-limited phase-preserving amplification can 
be implemented through non-degenerate parametric am¬ 
plification. In this process two distinct incoming modes, 
denoted here as the signal and the idler at the frequencies 
cjsi and cjid, are coupled to a strong pump mode by a 
nonlinear three-wave mixing element yielding amplified 
outgoing modes. We will first summarize the derivation of 
the amplifier input-output relations [21] before analyzing 
the erasure of the which-path information. 

A. Input-output relations 

We consider now a device with two ports and for each 
port we separate the incoming and outgoing modes. The 
device is operated in reflection but in the visualizations 
(see Figs. 1-2) we draw it in transmission for conceptual 
simplicity. For a reflective device, the signal input dsi,in(Q 
and output dsi,out(Q are related to the internal mode of 
the device dsi(t) through the coupling strength /^gi [22]: 

^si,out(^) ^si,in(^) — \/ ^si^si (^) • (1) 
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Because the internal mode is coupled to the external 
modes, it becomes damped at the rate and driven 
with V^^si,in(^): 

^si(^) ~ ^ [^si(^):-^a] \/^si^si,in(^)• (2) 

Similar equations hold for the idler mode aid(t). To solve 
for the outputs as a function of the inputs in Eqs. (l)-(2), 
we need to know the internal dynamics. For that pur¬ 
pose and for concreteness, we take the quantum-limited 
phase-preserving amplifier to be realized with a Joseph- 
son parametric converter amplifier (JPC) [23-25], whose 
Hamiltonian is a three-wave mixer, 

Ha/Ji= y] Wfcapfe-iAa^jatjap+iA*asiaida|^. (3) 

A:=p,id,si 

To operate the device as an amplifier, the pump mode 
Up is driven strongly at the frequency Q = cjgi + (^id such 
that it reaches a steady state ap(t) = o^pC^^^ The 

pump mode provides the energy for the amplification. By 
ignoring the remaining quantum fiuctuations in ap(t) and 
going into the frame rotating at the eigenfrequencies of 
the signal and the idler modes, the resulting dynamics is 
set by the Hamiltonian 

Hi/h= -iAa^ja/j + WagiaM, (4) 

where A = |A|e^^ = apX is the effective amplification 
strength. Using this Hamiltonian in Eq (2), the input- 
output relations (1) can be expressed in the frequency 
domain as (see Ref. [21] for details): 

^si,out(^) ~ ^(^)^si,in(^) T in(^)’ 

(5a) 

fiid,out(^) — 5^(^)d^id,in(^) T 6'^\/|^(^)P in(^)' 

(5b) 


The gain factor g{uj) and the amplification bandwidth D 
are 


9{^) = 
D = 


VG D ^ ^ /t D2 

1 ^id^si 

y/G T 1 ^id T ^si 


(6a) 

(6b) 



Figure 1. (a) Schematic of a quantum-limited phase-preserving 
amplifier (large triangle) and subsequent heterodyne detection 
of both quadratures; (b) An equivalent, high-gain representa¬ 
tion using a 50-50 beam-splitter (BS) and two phase-sensitive 
amplifiers (triangles). 


expressed in the frame rotating at the resonance frequen¬ 
cies; see Fig. 1(a). For simplicity we have ignored the 
phase between the idler and the signal ports. The next- 
to-leading order contribution is to ignore the second order 
terms but keep the first order terms in Eq. (6a), 

resulting in temporally non-local but causal responses 
with a delay kernel of the type e“^'^asi/id,in(^ “ r)dr 
in Eqs. (7). 


B. Erasure of the which-path information 

We see from the input-output relations (7), that when 
the idler input is the vacuum, both quadratures of the 
signal input are amplified along with an added extra half 
a quantum of noise originating from the idler. When 
both input ports contain signals, they are coherently 
superposed in the outputs and the which-path information 
of the signals is erased in the frequency domain. 

To further analyze the erasure of the which-path infor¬ 
mation, we specify that the signal output is measured by 
heterodyne detection [26, 27]: equal sampling of quadra¬ 
tures / (X Re(asi,out) and Q oc Im(asi,out), see Fig. 1(a). 
For simplicity we scale the signal with 1 / \/G and consider 
the high gain limit G ^ 1. The outcomes of a weak 
continuous measurement of infinitesimal duration dt are. 


where the total coupling rate k = Kgi -h /^id, the coupling 
asymmetry d,^ = (/^id — K,si)/i^ and the power gain \/G = 
15(0)1 = (fCsiKid + 4|A|^)/(KsiKid - 4|A|^). 

In the limit of a wide amplification bandwidth H, z.e., 
when the input signals are slowly changing with respect to 
the time scale the gain factor g{uj) (6a) is a frequency 
independent constant to leading order (the zeroth order) 
in (jo/D. This implies time-local input-output relations: 

asi,out(i) = VGasi,in{t) + VG- la|d,in(*)> 

Sid,out(i) = \/Gaid,in(i) + VG- lali,in(*)> (7b) 


dlr{t) = V2Re (Ssi.in + Sid,in) dt + dWi{t), (8a) 

dQr(i) = V2lm (asi,in - Sid,in) dt + dWQ{t). (8b) 

They consist of two parts: the part expected based on the 
prior knowledge of the system and the unexpected part 
(the ‘innovation’) dWpQ. Alternatively, dWpQ represent 
the quantum noise of the channels and are modeled by 
independent Wiener processes with variance dt [26, 27]. 

From Eq. (8), we see that the amplification-detection 
scheme is equivalent to a 50-50 beam-splitter followed 
by phase-sensitive amplifiers in both of the output arms 
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implementing two single-quadrature (homodyne) mea¬ 
surements in the / and Q directions, see Fig. 1. This 
interpretation also illustrates the nature of which-path 
erasure: observers of the output ports cannot know where 
the signals came from. In the high gain limit, the idler 
output is fully entangled with the signal output [28] con¬ 
taining no extra information. This can be understood as 
two-mode squeezing by the amplifier having effectively 
erased two of the four incoming quadratures. 

The input-output relations (7) of the amplifier are ex¬ 
pressed with the Hermitian conjugated input operators, 
a|d in ^si in • uoticc that they can be interpreted as 
complex conjugation in Eq. (8) when the amplifier is fol¬ 
lowed by the quadrature measurements. This is consistent 
with the physical picture of unidirectional information 
flow from the signal sources. The observer of the unidi- 
rectionally traveling signals can only make measurements 
whose measurement backaction to the system is expressed 
by the operators a^, that is, e.^., observations of discrete 
photon emissions or continuous leaking of the cavity field. 
The interpretation of the amplification-detection stage 
through beam-splitters and quadrature measurements 
gives practical means to handle components of cascaded 
quantum network with Bogoliubov transformations in 
their input-output relations [29-31]. 


III. CONCURRENT GENERATION OF 
REMOTE ENTANGLEMENT 


We now study a phase-preserving amplifier as a which- 
path information eraser to concurrently generate re¬ 
mote entanglement. The considered configuration con¬ 
sists of two transmon qubits inside separate, remote 
superconducting cavities, see Fig. 2. The cavities are 
driven through weakly-coupled input ports and monitored 
through separate strongly-coupled transmission lines that 
form the signal and idler ports of a Josephson parametric 
converter amplifier. In the frame rotating at the cavity 
driving = ujj — Aj and the qubit frequencies the 
Hamiltonian for a dispersively coupled qubit-cavity pair 
is 


= (Aj + a]aj + ej{t)a] + (9) 

where Xj denotes the dispersive coupling strength. When 
the cavities are driven at their resonance frequencies Aj = 
0, they build up symmetric qubit-state dependent phase 
responses. In the ideal case these responses are identical 
and when they are amplified with a high-gain quantum- 
limited phase-preserving amplifier according to Eq. (7) 
there is no which-qubit information left in the outgoing 
modes. This allows pure joint measurements of the signal 
sources, here the transmon qubits, with measurement 
back-action that projects to an entangled subspace. 



Figure 2. Schematic of a concurrent remote entangling experi¬ 
ment showing two dispersively coupled qubits inside separate, 
driven and monitored cavities. The left I-Q planes show the 
qubit state dependent coherent states. The right I-Q plane 
shows the structure of the signal output, Eq. (7), where no 
which-qubit information is available leading to the joint mea¬ 
surement . 


A. Stochastic master equation of the joint 
measurement process 


To analyze the configuration and the entanglement gen¬ 
eration in detail, it is modeled with a stochastic master 
equation [26, 27] (SME) for the density matrix p of both 
qubit-cavity pairs. We derive it by using input-output the¬ 
ory for cascaded quantum systems [22, 29, 30] (essentially 
^si,in = a/^^ 1 ^id,in = and representing the 

amplification-detection stage as a beam-splitter followed 
by two quadrature measurements as in Eq. (8). See de¬ 
tails of the derivation in Appendix A. The resulting Ito 
stochastic master equation for the monitored qubit-cavity 
pairs becomes 



dt 


+ f2\-D[VfJai)p+lv 





dt 


1 

71 

1 

71 


+ ^/r]2K2a2) P<^Wi 
H {i^rnKiCLi - i^r]2H2a2) /SdlTg, 


( 10 ) 


where dp = p{t-\-dt)— p{t). The first two rows describe the 
open quantum system dynamics. The dissipator terms 
V (c) p = cpc^ — ^ {c't'c, p} model the coupling of the 
cavity fields dj to the transmission lines with loss rates 
K,j . The dissipators for the qubits describe relaxation and 
pure dephasing with the rates F-^ = IjTl and 

The last two rows describe the measurement backac¬ 
tion, which updates the best estimate of the quantum 
state based on the new information dWpQ in the het¬ 
erodyne measurement of the signal output. It is repre¬ 
sented by innovation terms [27] %{£) p = [c — Tr (pc)] p -1- 
p [c^ — Tr (pc^)] that are linear in the measurement op- 
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erators c but nonlinear in the density matrix p. The effi¬ 
ciencies appearing in the innovation terms, describe 
the fraction of the information measured by the observer. 
The remaining fraction of the information is lost to the 
environments and averaging over it leads to dephasing. 
Here the efficiencies are pi = pfji and p 2 = VgVV 2 - They 
consist of the transmission coefficients fjj and the mea¬ 
surement efficiency p of the amplification-readout chain. 
With a finite amplification gain G, there is an asymmetry 
between the idler and signal outputs in Eq. (7). When 
measuring only the signal output, the associated loss of in¬ 
formation is represented with the inefficiency l — Pg = 1/G 
(in practice, G ^ 10 — 20 dB [1, 25]). 


B. Stochastic master equation for the monitored 
qubits 


Given that we consider driven and damped disper- 
sively coupled qubit-cavity systems, the cavity states 
can be assumed as a superpositions of the qubit depen¬ 
dent coherent states [32] {|<ai) |gj), |<a^) |ej)}, see Fig. 2. 
The coordinates follow the classical equations of motion, 
«ig(i) = W -i (Ai ± y - iy) see Fig. 3(c)- 

(d). With modern superconducting technology, one can 
achieve typical values of Ti ^ 100 /is [33, 34] for the 
decay times of the qubit. In this regime, the probability 
for a Ti relaxation event in either of the qubits can be 
assumed negligible small during the measurement time 
Tm ^ 1 /is. Then the qubits’ dynamics can be seen 

as frozen out, except for the measurement backaction. In 
this case, the cavity dynamics can be integrated out from 
the full SME (10) by using these time-dependent coher¬ 
ent states as an Ansatz [32, 35]. Reduction of the full 
SME (10) is very helpful since we are primarily interested 
in the qubits and the full SME is rather inconvenient to 
deal with due to its large Hilbert space. 

We approach here the integration of the cavity dynam¬ 
ics from a simple point of view. First, by exploiting the 
linearity of the innovation terms with respect to the mea¬ 
surement operators, we notice that the qubit-cavity pairs 
in Eq. (10) can be organized so that the two blocks look 
superficially decorrelated, except for the same stochas¬ 
tic ‘driving fields’ dWpQ. In practice, Eq. (10) encapsu¬ 
lates a technical source of correlation through the terms 
(dj) pdWi^Q and (dt) pdWpQ of the innovation operators. 
However, as these terms stem from the normalization 
condition, we can first uncorrelate the two blocks by sacri¬ 
ficing the normalization. Then, we make use of the single 
qubit-cavity results of Ref. [32] for integrating out the 
cavity individually on both qubit-cavity pairs, and finally 
add the normalization. We have verified these heuristic 
arguments with a numerical comparison and with an an¬ 
alytic calculation where we explicitly integrate out the 
cavity dynamics from the full SME (10) using a positive 
P representation [35, 36]. 

The resulting SME for the qubits’ density matrix pq = 


Trci 2P is 



+ + S2{t)al) pqdWi 

+ (i5i(i)(T] - iS2{t)al) PqdWQ, (11) 


where Oj(t) = XjRe is the ac Stark effect 

induced by the photons in the cavities. The photons 
leaking out from the cavity carry information about the 
qubits’ population [35] causing measurement induced de¬ 
phasing at the rates r;^(t) = Xjim {(ai(t) [(a|(t)]*}. The 
information is encoded into the distinguishability of the 
pointer states see Fig. 2-3. Therefore we define the 

complex measurement amplitude Sj {t) of the operator 
as. 


Sj{t) = - a|(i)], (12) 


whose real and imaginary parts are related to the mea¬ 
surements in the / and Q directions, respectively. The 
measurement rate is r^(t) = \Sj{t)\‘^. 

The outcomes of a weak continuous measurement of 
infinitesimal duration of both / and Q quadratures (het¬ 
erodyne) of the amplified signal, from Eq. (8), can be 
expressed as 


d/r(t) (ReiSi(t)d] H- ReiS2(t)d^) dt -f dWi 

v2 

H— ^KeUi{t)dt-\ — ^KeU 2 {t)dt, (13a) 

V 2 V 2 

dQr{t) {lmSi{t)al - lmS2{t)d-l) dt + dWq 

v2 

-1- -^lmUi{t)dt - ^lmU2{t)dt, (13b) 

V 2 V 2 


where the expectation values, for example (Re5i(t)d^), 
are taken for the instantaneous qubits’ density ma¬ 
trix Pq(t). Additionally, the terms involving Uj{t) = 
yJPjtZj [ai (t) + Q^i(t)] are not informationally meaningful 
since they only deterministically offset the signals, and 
will be ignored in what follows. Importantly, SME (11) 
and the measurement outcomes (13) show that hetero¬ 
dyne measurement of the output of a quantum-limited 
phase-preserving amplifier implements a pair of two-qubit 
measurements corresponding to the operators 

Ai{t) = Re5i(t)d^ + Re«S 2 (t)d^, (14a) 

AQ{t) = lmSi{t)al - lmS2{t)dl. (14b) 


The complex phase and the magnitude of the measurement 
amplitudes Sj{t) are tunable in situ by the cavity driv¬ 
ing. Thus, the measurement operators can be changed 
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continuously from a simultaneous separate readout of 
Ai oc al and Aq oc into a joint entangling readout 
Ai (X al ± . 


C. Balanced driving for perfect erasure of the 
which-path information from dissimilar sources 


To utilize the readout for remote heralded entangling, 
we drive the cavities at resonance Aj = 0 implying 
Im Sj{t) = 0. For the most efficient entangling read¬ 
out, we would like make both measurement amplitudes 
equal S{t) = Si{t) = S 2 {t) throughout the measurement— 
including cavity transients and unequal cavity-qubit pa¬ 
rameters. In an entangling readout, one does not want 
to gather any single-qubit information. Remarkably, the 
matching of the measurement amplitudes can be achieved 
by simple engineering of the drive amplitudes. This result 
improves the flexibility and scalability of the concurrent 
remote entangling scheme. We get the balanced driving 
amplitude. 


dbal) 


{t) = 


X 2 V 2 K 2 


K2Si{t) - Si{t) 


I4 + X2 


Si{t) 


(15) 

as a result of solving the input of the cavity 2 for a given 
output S 2 {t) = 5i(t), visualized in Fig. 3. Cavity re¬ 
sponses in the Q direction are left unmatched because they 
carry no information about the qubits’ population, see 
Figs. 2-3. However, the noise dWq needs to be recorded 
because it encodes the stochastic relative phase shift be¬ 
tween the qubits due to the unequal photon shot noise in 
each cavity [1]. 

As the erasure of the which-qubit information in the am¬ 
plifier output is made perfect with the balanced driving, 
the subsequent heterodyne measurement realizes a mea¬ 
surement of the joint operator For the initial state 

Pi = I++) (++I, where dx\-\-) = |+), the measurement 
backaction projects the system into a heralded entangled 
state with the success probability Ps = 1/2. 


IV. THE QUANTUM FILTER 

In the quantum filtering, one assumes that the mea¬ 
surement records dlr{t) and dQr{t) are known from the 
initialization up to a time t and then one would like to 
know the best estimate of the state of the open quantum 
system conditioned on this particular measurement record 
and the initial condition. The stochastic master equation 
gives the incremental update dpq from pq{t) to Pq(t -1- dt) 
given the new information dlFpQ(t) in the measurement 
records and knowledge of the system Pq(t). Naturally, 
a way to obtain the quantum filter is to solve the SME 
numerically for each measurement trajectory individually. 
However, this may generally be a computationally ex¬ 
pensive and slow task. Thus, an analytic quantum filter 
would be much more appealing. 






0 0.4 0.8 -1 0 1 

t (/is) x/^Hea (MHz)^/^ 


Figure 3. Dynamics of two dissimilar cavities in the 
frame rotating at the cavity frequencies; (a) Driving am¬ 
plitudes: unbalanced ei = £2 (black), balanced (yel¬ 
low) and (green); (b) Measurement amplitudes Si (yel¬ 

low) and S 2 (green) where dashed lines refer to driving by 
the balanced amplitudes; Cavity dynamics for the unbal¬ 
anced (c) and balanced (d) driving amplitudes: (red), 

(blue), (light blue) and y^Oe (orange). The 

parameters are /^i/27r = 8 MHz, £m/27r = 1 MHz, and 
Kijxi = X‘ 2 1 = 1 ^ 1 1 1^2 — 2. The dynamics in the Q direction 
(Imo) are left unmatched as informationally insignificant. 


We now consider SME (11) with the balanced real 
measurement amplitudes S{t) = Si{t) = 52(t), 


dPq = E ^ Pq] dt + [d-J - al, pq]dWQ 

j = l V 

2 

+ E (\/fiai) Pq + ip + Pqj dt 


+ ^n{al+al)XdWu 


(16) 


where we have explicitly written the backaction of the 
Q-measurement in the form of stochastic phase rotation 
(the second term). The measurement currents are 

dl,{t) = ^ (al + gI) di + dWi, (17a) 
v2 


dQr{t)=dWQ. 


(17b) 


In the following we take into account only the increased 
dephasing rate of the qubits by the Ti relaxation processes, 
r 2 = ^ but ignore its effect on evolution of the 

qubits’ population. This is justified by long typical Ti 
times, Ti ^ 100 /is, with respect to typical measurement 
time Tm ^ 1 ps. 

To derive the quantum filter, one needs to apply the Ito 
rule for changing variables in stochastic calculus [26, 27]; 
If dX = u{t)dt -I- a{t)dW and F(X, t), then 


dF = 





a^{t) d‘^F\ 
2 5X2) 


dF 

dt + a{t)—dW. 


(18) 

Based on this Ito calculus, we have found the analytic 
solution for SME (11) that expresses the two-qubit state 
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Pq(t) conditioned on an actual stochastic measurement 
record [/m(^), Qm(^)] and the initial state pi = |++) (++|. 
The time-dependent full solutions for the most important 
two-qubit Bloch coordinates are (the rest are shown in 
Appendix B) 


{ZZ) (t) 
{ZI) (t) 
{XX) it) 

{XY) (t) = 


e cosh/m(t) — 1 

e-A(t) C 0 sh/m(t) + l’ 

sinhI^{t) 

e-A(t)coshJ„,{t) + l’ 

_/^rm_(0 I IziZZt 

,0 V2^s'^s^ vt 2 = 1 ^ 2^; 

e-A(t) cos0+(t) + cos [Qm{t) — 0- 


(19a) 

(19b) 

(19c) 

W] 


e-A(t) cosh/ni(t) -h 1 

_/rm _(0 I izi^Zt A (0 I X ^2 j \ 

g y2r,sKs^ vt 2 = 1 ^ 2^; 

e-A(t) 0+(t) + sin [Qm{t) — ©- 


(19d) 

(^)] 


e-A(t) cosh/m(t) + 1 

f YT^(t) I 1 —A(t) I plj-N 

{XI) (t) = 2 e“V^^ + ^^+^ 2 ^j 


(19e) 


e 2 ^ cosh cos 


Qmjt) 


-01 


e-A(t) cosh/m(t) + 1 


The ac Stark effect induced rotation angles are denoted 
with Qj{t) = /q Oj(r)dr and 0±(t) = 0i(t) ± 02(0- 

We have used the notation of {psi^s)~^ = Yl‘j=i{Vj'^j)~^ 
and defined the combined measurement efficiency pt = 
ViV^/ivi V 2 — V 1 V 2 ) that has an important role in the 
analysis of the entanglement fidelity. 

We define A(t) = rm(T)dr = |5'(r)pdr as the ap¬ 

parent total information content recorded by the observer 
up to time t. The stochastic measurement records /m(t) 
and Qm{t) are weighted integrals of the measurement 
outcomes (17), 


Imit) + iQmit) = V 2 / S{t) [d/r(r) + id(5r('r)]. (20) 

Jo 

The measurement amplitude S{t) is the correct relative 
weighting function—a matched filter—between the dif¬ 
ferent time instances of the measurement. The analytic 
solutions of Eqs. (19) and (Bl) have been verified by com¬ 
paring them to the numerical solution of the stochastic 
master equation with perfect overlap within the accuracy 
of the numerical methods. 

Quite surprisingly, even for a lossless system, the co¬ 
herence of a two-qubit state is reduced by the factor 
exp(—rm(t)/27^s/^s) [3, 37-39] during the measurement 
[see Fig. 4(b)], where 



( 21 ) 

This reduction is due to an accumulated temporal mis¬ 
match between the acts that cause the qubit dephasing 
and the measurement backaction, at the rates and 
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Figure 4. The panels (a)-(b) show a quantum trajectory 
given by the quantum filter (19) corresponding to a sam¬ 
ple measurement record: {XX) (t) (blue), {ZI) (t) (ma¬ 
genta), {ZZ){t) (red), purity Trpq(t) (gray) and concur¬ 
rence C(t) (green). The analytic and numeric (not shown) 
solutions overlap perfectly. The most probable quantum tra¬ 
jectory [41] from the initial pi to the final state pf is shown as 
a dashed line; see Appendix C. Noise in the /m measurement 
carries information about the qubits’ populations and thus 
the measurement backaction affects all two-qubit Bloch co¬ 
ordinates. Noise in the Qm measurement carries information 
only about the relative phase between the qubits, causing 
fluctuations, e.p., in the {XX) coordinate without affect¬ 
ing to the populations or concurrence, (c) The two-qubit 
state of Eq. (19) shown as as a function of the measurement 
outcomes (/m,Qm). The parameters are /^i/27r = 5 MHz, 
Ki/xi = 1, « 2 /ki = X2/X1 = 1-1, = q = 0 and 

Pj = 1. The cavity driving ei{t) is on from ts — 0.03 ps 
to te = 0.5 ps [vertical lines in (b)] with the amplitude 
£m = 1 ^ 1 /V^ and a temporal shape similar to Fig. 3(a) re¬ 
sulting in Am = A(Tm) ~ 37r. The cavity 2 is driven with the 
balanced driving amplitude of Eq. (15). 


F^(t)/2, respectively. Physically, each qubit is entangled 
with cavity photons and the entanglement is not removed 
by the measurement backaction until the photons have 
leaked out. At the end of a measurement t = Tm, when 
the cavities have been brought back to the vacuum and 
all the available information has been recorded, the pu¬ 
rity revives. Interestingly, during cavity transients the 
dephasing rate r;^(t) can have negative values, which 
implies revival of qubit coherence. This non-Markovian 
dynamics originates from coupling the qubits to Marko¬ 
vian reservoirs indirectly through the cavities. The so¬ 
lutions (19)-(20) generalize and go beyond the previous 
results [1, 3, 32, 37-39] by deriving the quantum filter for 
a concurrent entangling two-qubit readout and verifying 
the purity reduction directly from stochastic calculus. 

The solution (19)-(20) is an analytic quantum filter 
without need for stochastic numerical solutions that would 
be limited by time step approximations [3, 40]. The 
quantum filter can be interpreted in two ways. First, 
given an actual measurement record [Im{t)^Qm{t)] it 
draws the stochastic quantum trajectories of the two- 
qubit state, see Fig. 4(a)-(b). In another interpretation. 
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the quantum filter gives the two-qubit state as function of 
measurement outcomes (/m,Qm) representing the effect 
of the measurement backaction, visualized in Fig. 4(c). 
The measurement outcomes Qm are normally distributed 
with zero mean and variance = 2A(t). For a strong 
measurement, the distribution of the measurement out¬ 
comes Jm is a 1 : 2 : 1 mixture of three normal dis¬ 
tributions, each with variance = 2 A(t), centered at 
7m = — 2 A(t), 0, 2 A(t). This gives a definition of the 
measurement strength A/m/cr = ^j2K{t) as the distin- 
guishability of the parity subspaces. 


A. Concurrence 


To examine the fidelity of the entanglement, we calcu¬ 
late the concurrence C [42] from the quantum filter solu¬ 
tions [see Fig. 4(b)]. In the limit of strong measurement 
A/m/cr 1, the state has collapsed with high probability 
either to an entangled state with {ZZ) = —1 or to a 
product state with (ZZ) = 1. In this limit, the concur¬ 
rence can be accurately approximated from the simplified 
expression [43]. ^^(t) — 2max‘(0, |pge,eg| 'y/Pgg^ggPee(eej"* 
Expressing this with the quantum filter solutions results 
in 


C{t) = max < 0 , 


exp 


( 377t-l A(t) 
V 2 


- Tit) - 1 


cosh Jm (t) -f exp A(t) 


( 22 ) 

where FI = F^ + F^ is the sum of qubits’ dephasing rates. 


Let us now consider the case of where the measurement 
has ended such that T^{t = T^) = 0 . Curiously, there 
exists an important bound for the total measurement 
efficiency r]t since the numerator of Eq. (22) needs to 
be positive for an entangled state with concurrence C > 
0. Even for ideal qubits (F^ = 0), there is a threshold 
77 t > 1/3 for forming a entangled state. For symmetric 
transmission this corresponds to 771 = 772 > 1/2. In the 
presence of decoherence the bound naturally becomes 
stricter: 


Vt > 


1 


Q O 


(23) 


Above the threshold, the purification by the measurement 
backaction dominates over the measurement induced and 
qubits’ natural dephasing. In addition one can see from 
Eq. (22) that given a fixed measurement time Tm and 
non-ideal efficiencies Tjj < 1 , there exists an optimal mea¬ 
surement strength Aq that maximizes the concurrence. 
Physically this can be understood as follows: with the op¬ 
timal Aq, the purification by the measurement backaction 
is in balance with the measurement induced and qubits’ 
natural dephasing. 


V. DISCUSSION AND CONCLUSIONS 

Compared to corresponding photon-counting based con¬ 
current remote entangling schemes [12, 13] whose entan¬ 
glement fidelity is more robust to losses and inefficien¬ 
cies, the proposed continuous variable scheme achieves 
very high generation rate ^ 10 ^ s~^ of entangled qubit 
pairs (Trep r\j 5 /is). Experimental values (l/F^ = 15 /is, 
Tm = 1 yus, r]j = 0.7) reachable in near-future supercon¬ 
ducting circuit experiments result in concurrence C ^ 20 %. 
The current experimental capabilities are such that the 
cavity-qubit asymmetries can be reduced through our 
pulse engineering scheme to the point that they will not 
be a limiting factor—rather it is the overall efficiency of 
transmission and measurement, where future technical 
improvements will lead to considerably better concurren¬ 
cies. 

In conclusion, we considered a readout chain of a 
quantum-limited phase-preserving amplifier followed by 
heterodyne detection and developed a physically intu¬ 
itive description compatible with the theory of cascaded 
quantum systems. Based on a stochastic master equa¬ 
tion approach, we theoretically demonstrated that the 
amplifier can be utilized as an eraser for the which-qubit 
information, even from dissimilar sources, and an element 
in a promising protocol for concurrent entanglement gen¬ 
eration between remote superconducting qubits. This 
protocol is feasible with existing technologies and can be 
expected to demonstrate formation of entangled remote 
qubits, primitive constituents of quantum communication 
and distributed quantum computation. 
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Appendix A: Cascaded quantum systems with a 
quantum-limited phase-preserving amplifier and 
heterodyne detection 

To rigorously derive the stochastic master equation 
of a cascaded quantum system [29], an effective method 
is to construct the unidirectional quantuni network by 
using the input-output triplets G = (S', T, H) of the net¬ 
work elements [30]. The G-triplet contains the scattering 
matrix S for the input-output ports of the element, the 
vector L that specifies the coupling to input-output ports 
and the internal Hamiltonian H. Eor compiling a net¬ 
work, one needs to know the rules for the cascade (series) 
G 2 < Gi and concatenation (parallel) Gi ffl G 2 prod¬ 
ucts. Let us consider two systems Gi = (Si,Ti,^i) and 
G 2 = (S 2 , £ 1 , 7 / 2 ), then the cascade and concatenation 









(Ala) 

(Alb) 


products are, respectively, 

G 2 <1 G\ = (S' 25 ' 1 , S 2 L 1 + L 2 , 

'Ll 


Gi ffl G 2 = 


Si 0 

0 S 2 


, Hi + . 


In the cascade product, the outputs of the system Gi are 
connected to inputs of the system G 2 , and the cascaded 
Hamiltonian has a corresponding driving term = 
H1+H2- '^{L\S 2 Li - l\sIL 2 ). 

The unidirectional quantum network of the concurrent 
remote entanglement setup of Figs. 1 and 2 is shown in 
Fig. 5. The G-triplets of the individual elements are: 


Gi 


Gi 


Gi 


G 




Gbs 



vT’ 

0 

0 _ 1 

~\/l ~ Vg 0 

/ 1 /V 2 0 0 
0 1 0 
0 0 1 
1 /V 2 0 0 


V- 



0 , 0 


0 , 0 


0 , 0 


(A2a) 

(A2b) 

(A2c) 

(A2d) 

(A2e) 


where G^ denotes the semi-classical driving of the j :th cav¬ 
ity, represents the internal dynamics of the j :th qubit- 
cavity pair with the Hamiltonian Hj = h{Aj + ^aDa^aj^ 

G^ models the losses of the transmission line, G^^^ de¬ 
notes the effective losses of the amplification-detection 
stage with r]g = {G — 1)/G, and Gbs denotes the effec¬ 
tive beam-splitter of the amplification-detection stage 
shown in Fig. 1. We write the identity (pass-through) 
element of a single G/ = (1,0,0) and several parallel 
ports Gj^^ = (/n,0,0). The non-monitored channels are 
denoted with a termination by a gray box in Fig. 5. 

Compilation of the network elements, in accordance 
with the subparts denoted with dashed boxes in Fig. 5, 
results in 


G^ = (G^ ffl Gi) < {Gi ffl G^) < (g^^ ffl Gj) , (A3a) 

G| = {Gi ffl G^ ffl Gi) < (g^ ffl Gf ^) < (g^ ffl Gf ^) , 

(A3b) 

Ga = (G, ffl Gbs ffl G,) < (g^^^ ffl G^,) , (A3c) 


where G^ denotes a driven qubit-cavity pair connected 
into lossy transmission lines and Ga is the effective 
amplification-detection stage by the quantum-limited 
phase-preserving amplifier and heterodyne detection of 
one of the output ports. The full compilation of these 
subparts gives the total network triplet 

G^ = Ga <1 {G^ ffl ~ Ht) • 





Gl 


Figure 5. The concurrent remote entanglement configuration 
of Fig. 2 in terms of the cascaded triplets G = (S', L, H). 
From left to right, the dashed boxes denote the compiled 
subparts G^, G^, and Ga- The non-monitored channels are 
denoted with a termination by a gray box. 


To form the stochastic master equation of the full 
quantum network, we need to know the Hamiltonian 
Ht: = J2jHj + hsj{t)a) + hs*{t)aj = J2jHj{t), which 
equals to Eq. (9), and the coupling vector 


Lt — 


y/Ki{l- T] i)ai \ 

^ {^/Kimai + ijK2mr]ga2) 


IF 


^ {i/nirtiai - 1 /K2ri2riga2) 




^K2{1 - ??2) a2 
i/k 2{1- Vg)a2 


J 


(AS) 


which have been simplified by absorbing the driving terms 
into the Hamiltonians. Finally, the stochastic master equa¬ 
tion corresponding to homodyne detection of the channels 
2 and 4, see Fig. 5, with the measurement efficiency 77 can 
be written as 


1 . ^ . 
dp= —[HT,p]<it + '^'D{Lk)pdt 


k=l 


+ U{i/vL2)pdWi + U{ii/vL4)pdWQ 

2 r ^ 

F ) ^^ (yT+^aj) p i di 


+ {^/piKiai + ^/ri2K2a2) pdWi 

+ T-H {iy/piKicii - ii/ri2K2a2) pdWq. (A 6 ) 
1/2 


In the main text, the small difference between Hj and 
tZj + is omitted for simplicity (or absorbed into the loss 
coefficients of the transmission lines fjj) and the qubits’ 
relaxation and pure dephasing terms have been added. 
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Appendix B: The quantum filter 


Here we list the rest of the two-qubit Bloch coordinates 
of the analytic quantum filter solution of Eq. (19) for 
SME (16): 

{YX){t)=e V2^s«s+ r,, 2 = 

Q-Ht) gi ]2 0^ (t) — sin [Qm (t) — 0- (t)] 


e-A(t) C0sh/m(t) + 1 

/ rm(t) I i-m A(t) -s ;^2 ^ \ 

(YY) (t) =e V2r,s«s^ vt 2 = 


(Bib) 


_P-A(t) 


cos 0+ (t) + cos [Qm {t) - 0- (i)] 


e-A(t) cosh/i„(i) + 1 

t rm(t) I A(t) I T^l+^ 

{XZ) (t) =2e 12.71-1+ .71 2 +^ 2 t) 


_A(t) j m 

e 2 sinh ^ ’ cos 


Qm(t) 


-01 


g-A(t) cosh Ji„{7) + 1 

f I 1~''72 A(t) I 7-12 + ^ 

(ZX) ( 7 ) =2e~y^^^~^ -+^2t) 


.A(t) / (J) 

9 cm h 




+ 02 


e-A(t) C0sh/m(t) + 1 

/ TmC*) I I-'.?! A(t) . ply.\ 

(yz)(t)=-2e r,, 2 +^ 2 ^) 


e 2 ^ sinh sin — 0^ 


e-A(t) cosh/m(t) + 1 

/ rm(t) I 1 — ''72 A(t) I p2.\ 

(ZY}(t)=2e V 2.72-2 + .72 2 2 ^) 


_Acq . , i^(t) . 
e 2 sinh ^ ’ sin 


Qmjt) 


+ 05 


e-A(t) cosh /m (7) + 1 


{YI) (t) 


_ f Tm (t) I 1 ~ •'71 A(t) _|_p2 . \ 

= _ 2e V^ryi^i VI 2 1 + 


e 2 cosh sin 


Qm(t) 

2 


01 


e-A(t) cosh/m(t) + 1 
{IX) (t) =2e V 2.72-2 + .72 2 + 2 ^) 

X — 


e 2 cosh cos 


%^+02 


e-A(t) cosh/m(t) + 1 

/ TmCt) I 1 — ''72 A(t) I p2.\ 

{IY){t)=2e V 2.72-2 + .72 2 2 ^) 


A(t) 

e 2 cosh 


X 


sin (^y^i + 02 


(Blc) 


(Bid) 


(Ble) 


(Blf) 


(Big) 


(Blh) 


(Bli) 


e-A(t) cosh/m(t) + 1 
The notation is the same as introduced in Sec. III. 


Appendix C: The most probable trajectory to an 
entangled state 

Given the initial pi = |++) (++| and final state 
p(^m) = Pf, as well as the time evolution in between 
described by SME (16), one may ask which are all the 


measurement trajectories {4, Qr} that connect the bound¬ 
ary conditions pi and pf, and what is the most probable of 
these trajectories, denoted here {/r,Qr}- First, the inver¬ 
sion of the quantum filter of Eqs. (19) and (Bl) gives the 
pair of values (/m, Qm) corresponding to the final state pf. 
Then Eq. (20) shows that the initial and final values are 
connected by all those trajectories {/r,Qr} that produce 
the pair (/m,Qm) via importance weighting (20). How¬ 
ever, this does not answer to the question what is the most 
probable trajectory. For that purpose, one needs to resort 
to more advanced methods in general. A possibility is to 
formulate a probability distribution of the measurement 
trajectories and apply the SME to make the connection 
between the time evolution of the quantum state and the 
measurement trajectories and, finally, by extremizing the 
distribution one finds the most likely trajectory [41]. 

Given the quantum filter of Eqs. (19) and (Bl), some of 
the most probable paths can be found by simple probabilis¬ 
tic arguments. Let us concentrate on the measurement 
in the I direction and consider the trajectories that 
end up in an entangled state for which (ZZ) = —1. In the 
strong measurement limit, the end points of the weighted 
trajectories /m(0 are normally distributed around three 
means —2A(t), 0, 2A(t). All the trajectories whose 

endpoint /m(Em) is around the origin, correspond to an 
entangled state (ZZ) = —1, see Fig. 4(c). Then the 
most probable end point producing an entangled state 
is Im{Tm) = 0- To see what is the most probable of all 
the trajectories Ir{t) that give rise to Im{Tm) = 0 we 
remember that the initial state is (ZI) (t = 0) = 0 and 
the weighting equation is, 

dl^{t) = V2S{t) {V2S{t) (ZI) {t)dt + dWi(7)) 

= V2S{t)dIr{t). (Cl) 

This shows that, in general, the probability cost of gener¬ 
ating an update of the quantum trajectory is minimized 
by choosing identically dlTi(t) = 0. For this update and 
the initial state |++) (++|, SME (16) keeps (ZI) (t) = 0 
and gradually updates (ZZ) towards (ZZ) = —1. Thus 
by taking this update dlTi(t) = 0 at every time step, the 
trajectory is /^(t) = 0 and it results in the most probable 
end point /m(Tm) = 0. Based on this, the most prob¬ 
able trajectory to the most probable entangled state is 
identically zero /^(t) = 0. 

For the measurement in the Q direction, the end point 
Qm(Tm) corresponding to the final state is solved again 
from the quantum filter of Eqs. (19) and (Bl). Then from 
the weighting equation, 

dQ„,(7) = V2Sit)dWQit) = V2Sit)dQrit), (C2) 

we see that the probability cost to reach the final value 
is minimized by updating the trajectory by dQr{t) = 
Q'^S{t)dt/(7‘^, where = Qm(Tm) mod 27r and = 
^/2/+52 (r)dr. This gives rise to the most probable 
trajectory Qr{t) = Q'^ 5'(r)dr/cr^. In Fig. 4(a)-(b), we 
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have shown the Bloch coordinates corresponding the most 


probable measurement trajectory as a dashed lines. 
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